Tunable Graphene Antennas for Selective Enhancement of THz-Emission 
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In this paper, we will introduce THz graphene antennas that strongly enhance the emission 
rate of quantum systems at specific frequencies. The tunability of these antennas can be used to 
selectively enhance individual spectral features. We will show as an example that any weak transition 
in the spectrum of coronene can become the dominant contribution. This selective and tunable 
enhancement establishes a new class of graphene-based THz devices, which will find applications in 
sensors, novel light sources, spectroscopy, and quantum communication devices. 
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I. INTRODUCTION 



The advance of optical antennas has led to incredi- 
ble new possibilities to control light-matter- interactions; 
including the directive emission of quantum systems or 
the modification of radiative rates at which such hybrid 
systems emit light i"— The parameter usually used to 
quantify modified emission rates of quantum systems is 
the Purcell factor F defined as F = 7rad/7rad with the 
radiative emission rate 7 ra d- The superscripts "a" and 
"fs" denote the presence of the antenna or an emission 
in free space, respectively. For a cavity, F is propor- 
tional to the quality factor Q and A 3 /V, where A is the 
wavelength and V is the mode volume. High Purcell 
factors have been achieved either in high-Q resonators, 
tolerating larger mode volumes, 4 or in systems support- 
ing localized surface plasmon polaritons (LSPP) with ex- 
tremely small mode volumes, therefore tolerating lower 
Q-factors* 5 - — However, since comparable Purcell factors 
appear for modes at higher frequencies in most of such 
implementations, these approaches are inappropriate for 
the selective enhancement of emissions at individual fre- 
quencies exclusively. Moreover, the impossibility to tune 
the spectral response of metal-based antennas over ex- 
tended spectral domains is extremely restrictive for many 
applications. 

Recently, graphene has been suggested as a possi- 
ble material that may remedy at least some of these 
limitations. 7 Moreover, since its operational domain 
would be in the THz frequency range (y « 0.1 — 30 THz), 
it may also be used for long-needed functional devices 
in the so-called THz gap. In this frequency domain, 
graphene offers low-loss SPPs^ — whose effective wave- 
lengths are two to three orders of magnitude smaller than 
the free space wavelengths at the same frequency^ In 
addition, the electronic properties of graphene can be 
widely tuned using the electric field effect . 13 ! 14 A final 
advantage of graphene is the much weaker interaction 
with light in the visible when compared to THz frequen- 
cies. This property is extremely useful for the selective 
enhancement of THz emission bands since optical emis- 




Figure 1: (color online) A graphene antenna to selectively 
enhance THz emissions: A molecule gets excited in the vis- 
ible/UV. Subsequentially, it strongly emits at a single fre- 
quency in the THz regime. The emission frequency can be 
tuned by adjusting the chemical potential \i c . 



sion bands are only marginally influenced by graphene 
antennas. 

A referential example for an application in the THz 
domain is the identification of specific molecules, e.g., in 
any transmission, fluorescence, or Raman spectroscopy 
experiments which rely on previously collected spec- 
troscopic information. 15 Although graphene-based op- 
tical antennas were already suggested to support such 
processes 7 -^, in principle, metal-based antennas can also 
be used in such applications. Therefore, the two key ques- 
tions to be answered are: How can graphene antennas be 
used to observe effects which are fundamentally inacces- 
sible by metal antennas and how can basic limitations of 
the latter ones be overcome? 

Here we introduce a concept that allows the selective 
enhancement emissions in molecules using tunable THz 
antennas made of graphene. It will be shown that dis- 
tinct and very weak spectral features at disparate fre- 
quencies can be enhanced without any structural mod- 
ification of the antenna. The scheme we are going to 
introduce provides the possibility to probe for all spec- 
tral features of a certain molecule; thus simplifying the 
spectroscopic classification of the molecule tremendously. 
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To enhance spectrally weak features, it is essential that 
the enhancement is selective, i.e. no other transitions 
shall be strongly enhanced except the desired one. This 
requires carefully designed antennas where higher-order 
resonances are shifted off the spectral domain of interest. 

For this reason, we will discuss the optical properties 
of a graphene-based antenna and elaborate on the details 
of the ultimate design first. All spectral properties are 
clearly motivated by the requirements described above. 
We will also discuss the interaction of the antenna with a 
specific molecule and outline how its spectral properties 
can be selectively enhanced. 



II. EFFICIENT AND TUNABLE GRAPHENE 
ANTENNAS FOR STRONG SPONTANEOUS 
EMISSION ENHANCEMENT 

Although graphene was suggested above as an ideal 
material, naturally its linear properties do not violate 
Kramers-Kronig relations and the dispersive properties 
are accompanied by absorption. Therefore, the rate at 
which radiation can escape an antenna-emitter system 
7 ra d is lower than the total emission rate 7 to t of the emit- 
ter, 7tot = 7ioss + 7rad- Both quantities are related by the 
antenna efficiency 77 = 7 ra d/7tot- The enhancement F to t 
of the total emission rate 7 to t of a dipole emitter in the 
vicinity of an antenna compared to the case of free space 
can be calculated within the Weisskopf-Wigner approxi- 
mation as F tot = 7t a ot/7t f ot = Ptot/Ptt - 17 Here, P t a f is 
the total power emitted by the dipole. Using the latter 
relation, it is possible to calculate the emission rates in 
the realm of classical electrodynamics. Assuming that 
the emission in free space is free of absorption (r/ s = 1), 
the Purcell factor with respect to the antenna reads as 
F = r] a F tot . In general, F to t depends on the position 
of the dipole, its polarization, and all properties of the 
antenna. 

In the following, we will show how to achieve a high 
Purcell factor for a graphene-based antenna in the THz 
regime at a preselected frequency exclusively. The elec- 
tric field effect can be used to change Graphene's chemi- 
cal potential \i c and consequentially its conductivity. See 
App. [A] for more information and details on our numer- 
ical implementation of the material. We further assume 
that the devices required to apply the electric field effect 
do not strongly influence the overall emission properties. 
This simplifies the analysis to free standing elements. 

The design of plasmonic antennas with high efficiency 
rf at a certain frequency strongly depends on geometrical 
parameters. Notably, the efficiency of sub- wavelength an- 
tennas can be quite low. To achieve a high Purcell factor, 
one thus encounters a trade-off between high F to t and r] a . 
Plasmonic antennas consisting of two elements generally 
exhibit a larger efficiency than single element antennas. 18 
In consequence, a configuration is considered where the 
dipole is situated in the center between the two elements 
of the antenna. The dipole is oriented in x-direction, 
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Figure 2: (color online) (a) & (c) Two circular graphene ele- 
ments cause a strong enhancement of the total emission rate 
Ftot for a dark antenna mode at A « 7.9 /mi (y « 38 THz). 
To visualize the field outside the antenna, the colormap has 
been truncated (white regions), (b) & (c) The antenna with 
two elliptical elements exhibits by far the strongest F to t for 
the desired dipolar mode at A ~ 12 /mi (y « 26.6 THz). 
The scale shown in (b) applies for all field plots with order-of 
magnitude rescalings. 



i.e. parallel to the connection line of both elements. As 
a result of the strong confinement of the LSPP, a small 
separation between the antenna elements is required. Ini- 
tially we chose an antenna design consisting of two circu- 
lar elements. Such a device can be fabricated e.g. using 
nanosphere lithography and the single discs may be de- 
scribed as Fabry- Perot resonators j 19 i 2Q A configuration 
with radii of 200 nm and separation of 30 nm was found 
to exhibit F to t ~ 10 6 -7 , and an efficiency 77 « 0.4 at 
v « 26.6 THz with \i c = 1 eV. However, this structure 
also permits the coupling of the dipole to higher-order 
resonances of the antenna, see Fig. [2] (a). These higher- 
order resonances, sometimes called dark modes, have ex- 
tremely low efficiencies since they are only weakly cou- 
pled to the far-field. Nevertheless, they exhibit a strong 
F t ot and are therefore undesired. For our present sce- 
nario, a dark mode at v « 38 THz had an enhancement 
exceeding that of the dipole resonance with an efficiency 
of only r] ~ 0.035. Therefore, these higher-order res- 
onances have to be suppressed in order to render the 
device useful. 

We found that the problem can be circumvented by 
stretching the discs in the vertical (y) direction, i.e. per- 
pendicular to the dipole, by a factor of two. This stretch- 
ing results in an increase of the effective radiating area of 
the antenna and thus of its efficiency. This is accompa- 
nied by a blueshift of the undesired higher-order modes. 18 
We observed a suppression of these dark modes by several 
orders of magnitude for an antenna made of two elliptical 
elements. Another aspect in favor of our design is a fur- 
ther enhancement of F to t, see Fig. [2](c). Generally, F to t 
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Figure 3: (color online) (a) A plane wave with amplitude Eo 
excites the dipolar resonance A ~ 12 fim of the antenna, (b) 
& (c): Distribution of |E (r) / Eo\ in the plane of the graphene 
ellipses (x-y-plane) and in-between at x = (y-z-plane). The 
region with strong enhancement is approximately 30 x 100 x 
20 nm 3 (blue region). For a dipole situated here, the Purcell 
effect is comparable to a placement in the center. 

at other frequencies than the antenna's dipole resonance 
is two or more orders of magnitude smaller. This offset 
permits the desired selective enhancement of molecular 
resonances if the resonance of the antenna coincides with 
one of the resonances of a quantum system as it will be 
outlined in the following section. 

One may furthermore ask the question how strongly lo- 
calized the Purcell enhancement is for the found efficient 
antenna. Because of the plasmonic nature of the mode, 
one expects a rather small region of enhancement. To 
measure the size of this region, one can use the reciprocity 
theorem : 21 ! 22 Illuminating the antenna with a plane wave 
corresponds to an emitting dipole at infinity. We assume 
that the plane wave excitation comes from a direction of 
strong emission of the antenna's dipole mode. Then, the 
field enhancement at some position is a measure for the 
radiated power of a dipole at that very position. We fur- 
ther assume that the dipole emitter couple to the dipole 
mode of the antenna. Then, radiation patterns do not 
change strongly for different emitter positions. Conse- 
quentially, the ratios of field enhancements are propor- 
tional to the ratios of the overall radiated power and thus 
the ratio of Purcell enhancements. We have made a cor- 
responding calculation. As expected, we found a region 
of strong interaction confined in-between the antenna, see 

Fig. m 

III. THE CHANGED SPECTRUM OF 
CORONENE 

Up to this point, the enhancement of the total emis- 
sion rate F to t was referring to a dipole in the vicinity of 
a graphene antenna. To demonstrate the performance 
of the suggested implementation, coronene (C24H12) has 
been chosen as an example. It belongs to the group of 
polycyclic aromatic hydrocarbon (PAH) complexes and 
is of both great practical and theoretical interest. 23 The 



molecule exhibits resonances in the THz regime and its 
properties are well-documented. 24,25 

The emission process of coronene can be understood as 
a stepwise procedure following Ref. [23|. First, upon ab- 
sorption of a UV photon, the electronic configuration of 
the molecule is transferred into an excited state. The ex- 
citation will be assumed in the near ultraviolet around a 
frequency of ~ 990 THz (300 nm) where coronene absorbs 
strongly. The absorbed energy is then almost instanta- 
neously redistributed to the different vibrational modes 
of the molecule. The characteristic time scale of this pro- 
cess is about 10 -12 s. This redistribution obeys a Planck 
distribution B (co>;T m ) over the respective eigenenergies 
hwi. The temperature T m refers to the microcanonical 
temperature of the molecule and can be approximated 
by T m ~ 2000 (E (eV) /N c f A K with the energy of the 
exciting radiation E in electronvolts and the number 
7V C = 24 of carbon atoms of coronene. The expression for 
T m can be derived from the density of states of PAH's. 
It has been proven correct within an error of ~ 10% for 
35 . . . 1000 K. 23 For coronene and the given excitation, we 
get T m « 990 K. The emission spectrum in free space is 
then given by the oscillator strengths / (u^) weighted by 

Close to the antenna, the emission of the molecule is 
strongly modified. We assume here, that the redistribu- 
tion of energy into the different eigenmodes is much faster 
than any involved radiation process. This seems reason- 
able since the natural emission time for THz radiation is 
in the order of or even bigger than one second. 23 There- 
fore, one might hypothesize that even an enhancement of 
several orders of magnitude for the given THz transitions 
in coronene should be possible without any noticeable in- 
fluence on the internal redistribution dynamics. Hence, 
by enhancing the emission rate for just one eigenmode, 
the emission of the remaining ones decreases. Then, one 
finds different total emission rate enhancements for each 
transition J^ot (^i)? which deviate from the Purcell fac- 
tors of a single dipole, i.e. J^ot (^i) 7^ ^tot (^)- For a 
linear multiresonant system, modified emission rates are 
related to each other by a sum rule. 26 Because of the in- 
ternal redistribution dynamics of coronene, this approach 
cannot be used. We assume that the power emitted upon 
excitation does not change if the antenna is present, i.e., 

P tot = T,i fi ^ni J7 ( UJ i) B ( UJ k]T m ) = P^ t implying en- 
ergy conservation. Then, the enhancement of the radia- 
tive rate with respect to the antenna efficiency is given 
by 

^rad (uJi) = T] (Ui) ■ F to t (^i) ' 

Efe ^fe7t f ot (^fc) B (cj fc ; Tm) ^ 
J2 k fi^kFtot (uk) 7tot (^fc) B(u k ]T m ) 

Since the characteristics of the graphene antenna can be 
tuned by changing the chemical potential /i c , the en- 
hancement also depends on fi c . If now a certain transi- 
tion is selectively enhanced to a large extent, F tot (ojj) ^> 
^tot (wk) (k ^ i), the energy supplied to the molecule 
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Figure 4: (color online) The modified emission spectrum of 
coronene in the vicinity of the graphene antenna as a function 
of the chemical potential /x c . For reference, the green curve 
displays the emission spectrum of the bare molecule exhibit- 
ing several peaks with different strengths. When tuning the 
chemical potential by taking advantage of the electric field 
effect, the antenna may strongly and selectively enhance only 
a single emission line, so that most of the energy leaves the 
molecule via this transition. 



IV. CONCLUSION 

In conclusion, we have suggested and verified that due 
to the unique properties of graphene, antennas made 
from this material can be used to selectively enhance dif- 
ferent molecular transitions at THz frequencies by several 
orders of magnitude. The achievable efficiencies and the 
suppression of undesired modes were found to strongly 
depend on the actual antenna geometry. A suitable im- 
plementation which consists of two closely placed ellip- 
tical graphene elements was introduced. Using this an- 
tenna, weak transitions can be made dominant in the 
emission spectrum. The potential of graphene antennas 
to selectively enhance molecular transitions can clearly 
be anticipated to lead to new applications in THz spec- 
troscopy but will also contribute to the development of 
novel sensors, highly directed single photon light sources, 
and quantum communication and quantum computation 
devices. 28 



will predominantly leave it by exactly this transition. 
This result would also hold with respect to the sum 
rule in Ref. [26|, for which we would find J 7 ^ 1 " (u)i) = 
Ftot (wi) I ^2 k Ftot (<Uk)- Hence, it might be expected that 
the shift of emission processes in favor of a single transi- 
tion by a selective enhancement is likely to be a generic 
property of molecular systems. The spectral resolution 
of the selective enhancement is limited by the antenna. 
It is naturally linked to the width of the antenna reso- 
nance which is approximately the relaxation rate T. So, 
if T is smaller than the minimum spectral distance be- 
tween certain resonances one wishes to distinguish, the 
selective enhancement should work. This is the case for 
coronene; current fabrication techniques allow relaxation 
rates in the order of a few meV or THz, respectively^ 

To demonstrate this selective and tunable enhance- 
ment, Fig. [H shows how the emission spectrum of 
coronene changes if the antenna resonances are tuned. 
Although extensive simulations were performed, only a 
few selective results are shown where the resonance was 
tuned to enhance the emission at frequencies correspond- 
ing to "jumping-jack"-, "drumhead"- and a "C-H-out-of- 
plane" -vibrational motion at 11.3 THz, 16.5 THz and 
26.0 THz, respectively. 27 The quantitative visualization 
of the spectra is performed by assuming, as usual, a 
certain line width for each resonance. For PAH's, we 
used a linewidth of 0.15 THz, which is characteristic for 
molecules isolated from each other. 27 It can be clearly 
seen that by tuning the antenna it is possible to choose 
the transition which is forced to emit light into the far- 
field; the main goal of this contribution. Various effects 
might broaden the lines as e.g. the coupling to surfaces 
or collisions with other molecules. Nevertheless, it can 
be safely anticipated that the general characteristics of 
the spectrum survive. 
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Appendix A: Numerical implementation of graphene 

At temperature T, the in-plane conductivity of 
graphene can be approximated using the Kubo formula 



a s (lj) i 
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with the chemical potential \i c and the relaxation rate 
Tjlii 2 - T was taken to be 0.1 meV, consistent with values 
reported earlier. 29 - Note that the second interband term 
is not temperature dependent for k B T <C |/i c |. This is in 
very good approximation for the values of the chemical 
potential used in our analysis, even at room tempera- 
ture, where k B T w 26meV. Throughout, we assumed a 
time dependency according to exp (— icut). In our calcula- 
tions, we represented graphene as a thin conductive layer 
with thickness d = 1 nm <C A and relative permittivity 
Eg (id) = 1 + ia s (uj) I (eoiod), see also Refs. M ll2ll30l - [32l 
and [29|. It was verified that none of the results are quan- 
titatively affected by that finite thickness. All simula- 
tions were performed using the Finite-Element Method 
as implemented in COMSOL. 33 
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